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ABSTRACT 

We study linear polarization due to scattering of light on a cloudlet of particles, taking 
into account the radiation drag and the gravitational pull exerted on them by a central 
body. Effects of special and general relativity are included by connecting a model of Bc- 
loborodov (1998) for the local polarization of scattered light with Abramowicz, Ellis & 
Lanza (1990) formalism for the particle motion near an ultra-compact star. Compact- 
ness of the central body and its luminosity are two critical parameters of the model. 
We discuss the polarization magnitude of photons, which are Thomson-scattered into 
direct and higher-order images. Importance of the latter is only moderate under typ- 
ical conditions, but they may give rise to distinct features, which we explore in terms 
of a toy model. The scattered signal exhibits variations of intensity and of polarization 
with mutual time-lags depending on the beaming/focusing effects and the light travel 
time. 
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1 INTRODUCTION 

Scattering of ambient light by fast moving flows has been 
recognized as a conceivable mechanism operating in differ- 
ent classes of objects. This process quite likely contributes to 
the linear polarization of initially unpolarized soft radiation 
up-scattered in blazar jets (Begelman & Sikora 1987), winds 
from accretion discs (Beloborodov 1998) and in gamma-ray 
bursts (Shaviv & Dar 1995; Lazzati et al. 2004; Levinson & 
Eichler 2004) . A significant level of intrinsic linear polariza- 
tion of IT ~ 3% was reported in a microquasar GRO J1655- 
40 (Scaltriti et al. 1997) and similarly for LS 5039 (Combi et 
al. 2004) in the optical band. A model of magnetized fireball 
polarization was discussed by Ghisellini & Lazzati (1999), 
who demonstrate that the expected polarization lightcurve 
should exhibit two peaks. Variations of the Galactic Cen- 
tre linear polarization were reported in the millimetre band 
(Bower et al. 2005). Recently, Viironen & Poutanen (2004) 
brought the attention to strong-gravity effects in polarime- 
try of accreting millisecond pulsars. 

The idea of X-ray polarization studies providing clues 
to the physics of accreting compact objects was discussed in 
seminal papers (Angel 1969; Bonometto et al. 1970; Light- 
man & Shapiro 1975; Rees 1976; Sunyaev & Titarchuk 1985). 
Here we study a simple model with Thomson scattering on 
electrons ejected outwards from the centre or falling back; 
a novel point is that we consider strong gravity effects. The 
source of seed photons can be identified with the surface of a 



central star lying in arbitrary (finite) distance from the scat- 
terer. Alternatively, it can represent an axially symmetric 
accretion disc near a black hole, a quasi-isotropic (ambient) 
radiation field, or a combination of all these possibilities. 
In order to avoid numerical complexities we do not discuss 
other radiation mechanisms that also produce polarization, 
namely, we do not consider synchrotron self-Compton pro- 
cess, which is the most likely process wherever magnetic 
fields interact with relativistic particles (see Poutanen 1994; 
Celotti & Matt 1994 for discussion and for further refer- 
ences) . 

In the non-relativistic regime, a conceptually similar 
problem was examined by Rudy (1978) and Fox (1994), 
who considered polarization of light of a finite-size star due 
to scattering on free electrons in a fully ionized circumstel- 
lar shell. These works showed non-negligible solid angle of 
the source subtended on the local sky of scattering parti- 
cles has a depolarizing influence on the observed fractional 
polarization. This result was discussed by various authors, 
because the scattering of stellar continuum is a probable 
source of net polarization of visible light in early-type stars, 
for which polarimetry has become a standard observational 
tool (e.g. Poeckert & Marlborough 1976; Brown & McLean 
1977). General relativity effects are negligible in these ob- 
jects, however, even the Newtonian limit is complex enough: 
the observed signal does not easily allow to disentangle the 
contribution of a circumstellar shell, the primary radiation of 
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the star and the effect of interstellar medium. In our case the 
situation complicates further, because the system is highly 
time-dependent and mutual time delays along different light 
rays must be taken into account. 

The general relativity signatures, which we discuss in 
this paper, reach maximum in a system containing a com- 
pact body with radius less than the photon circular orbit, i.e. 
R* < r p h- This can arise either as a result of a non-stationary 
system with R+^R+ty), presumably during a gravitational 
collapse, or it can represent a static ultra-compact star if 
such exist in nature. Here we impose spherical symmetry of 
the gravitational field, and so both situations are identical 
as far as the form of spacetime metric is concerned. The 
Schwarzschild vacuum solution describes the external grav- 
itational field of all types of compact objects within general 
relativity, provided that their rotation is negligible and self- 
gravity of accreting matter does not contribute significantly 
to the gravitational field. We will assume that these con- 
straints are fulfilled. A sequence of R+ — const situations 
can be employed in order to model a collapsing case. 

The formalism that we apply works for a system with 
arbitrary compactness, even if the case of i?* = const < r p h 
seems to be an unrealistically large compactness per se, 
likely violating the causality condition for neutron stars. 
Such a high compactness would thus be normally excluded, 
however, the situation has not yet been definitively settled 
(see Lattimer & Prakash 2004 for a recent review). Accord- 
ing to astrophysically realistic equations of state, neutron- 
star sizes do not reach ultra-compact dimensions; their typ- 
ical radii should exceed the photon circular orbit (Lattimer 
& Prakash 2001; Haensel 2003), and so the gravitational ef- 
fects are constrained accordingly. However, the possibility of 
Ri, being slightly less than r p h persists, and there seems to 
be a growing awareness now of the fact that the very high 
density regime needs to be explored further. 

The interest in ultra-compact stars has been recently 
revived mainly in connection with gravitational waves in 
general relativity (Chandrasekhar & Ferrari 1991; Kokko- 
tas et al. 2004) . One may consider also more exotic options. 
Strange quark stars (e.g. Alcock et al. 1986; Dey et al. 1998) 
can have their radii extending down to the Buchdal limit 
(see Weber 2005 for a recent review on other forms of quark 
matter and their relevance for compact stars). If nucleons 
can be confined at the density lower than nuclear matter 
density, then Q-stars could exist with a relatively high mass 
of ~ 1O 2 M0 and the radius as small as R+ ~ 0.9r p h (e.g. 
Miller et al. 1998). On a more speculative level is the idea of 
gravastars (Mazur & Mottola 2004), compactness of which 
can exceed the Buchdal limit. Options for detecting the ther- 
mal radiation from different kinds of ultra-compact stars 
have been recently discussed by McClintock et al. (2004). In 
the case of transiently accreting neutron stars, the crustal 
heating has been nominated as one of relevant mechanisms 
generating thermal emission from the surface (Haensel & 
Zdunik 1990; Brown et al. 1998). 

The issue of realistic equations of state for ultra- 
compact star matter is beyond the scope of the present pa- 
per, likewise the actual observational information on masses 
and radii of compact stars (see e.g. Haensel 2003). It will 
be convenient to introduce a dimensionless parameter, ( = 
1 — R*/r, which maps the whole range of radii above the 
star surface onto (0, 1) interval, so that the form of graphs 



does not change with the star compactness. We consider the 
whole range of ^ £ ^ 1 allowed by general relativity. It 
is worth noticing that the formalism used below could be 
readily applied also to the case of an accretion disc as a 
source of light near a black hole. Obviously, a black hole 
represents a body with the maximum compactness and, at 
the same time, it is the most conservative option for such 
an ultra-compact object (lower, axial symmetry of the disc 
radiation field limits the possibility of exploring the problem 
in an analytical way). 

Our paper takes general relativity effects into account, 
including the effect of higher-order images if they arise. Al- 
though the signal is usually weak in these images, favourable 
geometrical arrangements are possible and, even more im- 
portantly, photons of the higher-order image experience a 
characteristic delay with respect to photons following a di- 
rect course. This delay (examined in detail by Bozza & 
Mancini 2004; Cadez & Kostic 2005) could help revealing 
the presence of strong gravitational field in the system. 

Our model provides a useful test bed for astrophysically 
more realistic schemes. In the next section we formulate the 
model and describe calculations. Then we show comparisons 
with previous results of other authors. We build our discus- 
sion on the approach of Beloborodov (1998, for polarization) 
and Abramowicz et al. (1990, for the motion in combined 
gravitational and radiation fields in general relativity). In 
these papers the individual components of the whole picture 
were treated separately, while we connect them together in 
a consistent scheme. A reader interested only in the main re- 
sults on polarization of higher-order images from light scat- 
tered on a moving cloudlet can proceed directly to section 4. 

As a final point of the introduction, it is worth notic- 
ing that Ghisellini et al. (2004) propose a model of aborted 
jets in which colliding clouds and shells occur very near a 
black hole and are embedded in strong radiation field. Ac- 
cording to their scheme, most of energy dissipation should 
take place on the symmetry axis of an accretion disc. This 
would be another suitable geometry, in which a fraction of 
light is boosted in the direction to the photon circular or- 
bit and eventually redirected to the observer. One may fear 
that the Thomson scattering approximation is not adequate 
to describe a turbulent medium in which electrons become 
very hot (Poutanen 1994), however, the accuracy should be 
still sufficient for energies of several keV, at which planned 
polarimcters are supposed to operate. Also the seed pho- 
tons have different distribution when they originate from 
an accretion disc, but we examined the variation of the po- 
larization that an observer can expect from this kind of a 
system (Horak 2005), and our calculations confirm that the 
expected polarization has magnitude smaller but roughly 
similar to values predicted by a simple model adopted here. 

2 THE MODEL AND CALCULATIONS 

2.1 The set up of the model and reference frames 

We assume fully ionized optically thin medium distributed 
outside the source of seed photons. These primary photons 
follow null geodesies until they are intercepted by an elec- 
tron, which itself is moving under the mutual competition 
between gravity and radiation (we do not consider the ef- 
fect of magnetic fields on the particle motion and radiation 
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in this paper). Hence, we adopt the approximation of single 
scattering and we assume Schwarzschild metric for the grav- 
itational field of a compact body. We consider frequency- 
integrated quantities. Polarization vector of scattered light is 
propagated parallelly through gravitational field to a distant 
observer and, as consequence, the polarization magnitude 
n(r,n) and the redshifted intensity I = (1 + Z)~ 4 I(r, n) 
(expressed here in terms of the redshift Z) are invariant. 

Polarization is described in terms of Stokes parameters 
I, Q, U and V (Chandrasekhar 1960; Rybicki & Lightman 
1979): I has a meaning of intensity along a light ray, Q and U 
characterize the linear polarization in two orthogonal direc- 
tions (say ex and ey) in the plane perpendicular to the ray, 
and V is the circularity parameter. The polarization angle is 
defined by tan 2\ = U/Q. It gives the orientation of major 
axis of the polarization ellipse with respect to e x ■ One can 
form a polarization basis in the local space by supplement- 
ing ex and ey with another unit vector, ez, pointed in the 
direction of the light ray. Three parameters are necessary 
to describe a monochromatic beam, for which the condition 
I 2 — Q 2 + U 2 + V 2 holds. In case of partially polarized light 
the whole set of four parameters is generally required. It is 
then customary to define the degree of elliptical polarization, 
n = I' 1 (Q 2 + U 2 + V 2 ) 1 ' 2 , which satisfies < II < 1. 

A four- dimensional tetrad can be constructed by ex- 
tending the three base vectors and supplementing them with 
a purely timelike four- vector. A suitable choice of the tetrad 
is described below. Let us consider a simple case when the 
incident radiation field is axially symmetric in the labora- 
tory frame, LF, with the basis (fit, e x , e y , e z )} Indices of 
four-vectors with respect to a local-frame basis are manipu- 
lated by the flat-space metric, diag(— 1, 1, 1, 1). 

We orient e z along the symmetry axis; other two spa- 
tial vectors, e x and e y , lie in a plane perpendicular to e z . 
Further, we assume that scattering electrons are streaming 
along the symmetry axis with four- velocity u — v} e t + u z e z , 
where u = cy and u z — cy/3 (7 is Lorentz factor, /3 is veloc- 
ity in LF divided by speed of light). Later on we carry out a 
Lorentz boost to the co-moving frame (CF) of the scatterer, 
(et,e x ,e y ,e z ), which is equipped with timelike four-vector 
e~t — u and three space- like four- vectors e x — e x , e y — e y . 
Spatial part of e z is oriented in the direction of relative ve- 
locity of both frames. 

Each incident photon of the ambient unpolarized radia- 
tion gets highly polarized when scattered by a relativistically 
moving electron. The total polarization is eventually ob- 
tained by integrating over incident directions and the distri- 
bution of scattering electrons. In order to describe propaga- 
tion of scattered photons, we denote four- vectors n = p/p* 
(with respect to LF) and h = p/p* (with respect to CF), 

1 Hereafter we understand three-vectors as spatial projections of 
their corresponding four-vectors (and we do not introduce special 
notation for them; there is no danger of confusion). We stick to the 
conventional formalism of Stokes parameters, but we remark that 
it can be recast by employing a covariant definition of the polar- 
ization tensor, components of which are assembled using suitable 
combinations of Stokes parameters (Born & Wolf 1964). It has 
been argued (Portsmouth & Bertschinger 2005) that the latter 
approach may be found more elegant and useful for discussing 
the radiation transfer of polarized light through the medium in 
general relativity. 



where p is the photon four- momentum (a null four- vector) . 
Due to the axial symmetry we can assume n v = n y = 0. 

In addition to the above-defined reference frames LF 
and CF, we introduce two 'polarization' frames: the lab 
polarization frame (LPF) with basis (fit, ex, ey , ez), and 
the co-moving polarization frame (CPF) with the basis 
(§t, e x , ey, ez). LPF is defined in such a way that ez is 
the three-space projection of the propagation four-vector n, 
ex lies in the (e x , e z )-plane, and ey is identical with the LF 
tetrad vector e y . CPF is defined analogically and indicated 
by bars over variables. Our definition of the reference frames 
is apparent from figure 1. 

2.2 Stokes parameters in terms of the incident 
radiation stress-energy tensor 

We start by calculating the polarization of the scattered ra- 
diation in CPF. Conceptually the model of local polarization 
is equivalent to the one employed by Beloborodov (1998). 
The incident radiation is unpolarized with intensity h. It 
can be imagined as a superposition of two parallel beams of 
identical intensities, l[ = — T\j2, propagating along 
■hi four- vector. The two beams are completely linearly polar- 
ized in mutually perpendicular directions and the scattered 
radiation is a mixture of both components. In the adopted 
choice of reference frames (see the figure 1), the spatial pro- 
jection of the propagation vector n is identical with the spa- 
tial projection of ez- Unequal contributions J^ 1 ' and J^ 2 ' to 
the total intensity I result in net linear polarization of the 
scattered beam. According to this, I, Q and U are non-zero, 
whereas the circularity parameter V vanishes. 

Assuming that each scattered photon experiences one 
scattering event in an optically thin medium (r <JC 1), non- 
zero contributions to Stokes parameters are (Chandrasekhar 



1960) 

SI = AFi(l + cos 2 a;) , (1) 

SQ = — AIi cos 2^5 sin 2 w, (2) 

SU = — AIi sin2<p sin 2 u>, (3) 



where ui is the scattering angle between n\ and n, A = 
3r/(167r). The scattering takes place in the plane that forms 
an angle ip with x-axis. Angles ip and to can be expressed 
using direction cosines, which are defined here as spatial 
components of the propagation four-vector fn of the incident 
beam, i.e. nf — cos ip sin uj, n( = sin p sin u> and nf = 
cosw. We obtain 

SI = A(l +nfnf^ 7i, (4) 
8Q = A (jiYnY ~ nfnf^j 7), (5) 
SU = ~2An?nYli. (6) 

This form is useful, as it allows integrating conveniently the 
partial contributions over incident directions to obtain 

I = Ac (f u + f zz ^ , (7) 

Q = Ac (f YY - T xx ) , (8) 

U = ~2AcT XY , (9) 

for the total Stokes parameters of scattered light. We de- 
noted 
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Figure 1. Geometry of the problem and the definition of angles. 



rpfl 



n^n" Ii(rii) dQ 



(10) 



the stress-energy tensor of the incident radiation field. 

We remind that the incident radiation was assumed 
axially symmetric in the CF, therefore the only non-zero 
components in this frame are T u , f tz , f zz , f xx , &nd f yy . 
These are further constrained by symmetry: T xx = T vv — 
(T" - T zz )/2. A relation to CPF components can be ob- 
tained by rotation about j/-axis by angle Using equations 
(7)-(9) we find 

I = ±Ac[ (3T" - f zz ) - (f u - 3f zz ) cos 2 #] , (11) 

Q = \Ac (T" — 3T ZZ ) sin 2 #. (12) 

The Stokes parameter U vanishes due to axial symmetry. 
The scattered radiation is partially polarized either in the 
(x, z)-plane, or perpendicularly to it. The former and the 
latter case will be referred as longitudinal and transverse 
polarization, respectively. Later on, in Sec. 3, we will demon- 
strate that this change can occur also with a group of cold 
electrons, whose bulk motion is determined by the radiation 
and gravitational fields of a star and recorded in the lab 
frame. 

The degree of polarization is calculated directly from 
definition: 

m = U = in^lsin 2 ^ _ m 



i-n n 



where 



n m = 



3T Z 



3T* 



(14) 



This result is equivalent to eqs. (4)-(5) in Beloborodov 
(1998), who applied the model of Thomson scattering 
to winds outflowing from a plane-parallel disc slab. Be- 
loborodov also calculated the polarization of scattered ra- 
diation and he found that its sign depends on the wind ve- 
locity. In our notation the change is captured in IT m , which 
acquires values in the range — 1 ^ n m ^ 1. Its meaning is 
evident from eq. (13): the absolute value |II m | is the max- 
imum degree of polarization of the scattered light and the 
sign of II m determines the sign of Q-P arame ter. 

In order to determine the polarization magnitude as 



seen by an observer in LF we carry out the Lorentz boost 
(e.g., Cocke & Holm 1972; Rybicki & Lightman 1979). The 
angle of observation 3 is transformed according to: 



■ T> sin i9, cos $ = yT> (cos # — 0), 



(15) 



where T> = 7 _1 (1 — /3 cos i9) _1 is the Doppler factor. Stokes 
parameters are transformed in the same manner as the ra- 
diation intensity and the boost retains the four-vector e y 
unchanged: / = V 4 I and Q = V 4 Q. It follows that the 
polarization magnitude |II m | is Lorentz invariant. By trans- 
forming all relevant quantities to LF, we obtain Stokes pa- 
rameters of scattered radiation, 



\AcV 6 



I = 



7 2 [(1 - 3/3 2 )T" 
+4/3T tz - (3 - /3 2 )T ZZ ] sin 2 ■&, (16) 
AcD^ 2 [(1 + f3 2 ) (T" + T zz ) - 4/3T iz ] + Q. (17) 



2.3 Critical velocities 

The aim of this subsection is to connect, in a self-consistent 
manner, the properties of particle motion through the ambi- 
ent radiation field with Stokes parameters of scattered light. 
In order to prepare for this discussion it is useful to intro- 
duce two critical velocities of the particle motion. 

Firstly, of particular interest is the velocity at which the 
polarization of scattered radiation vanishes (Beloborodov 
1998). The condition for velocity follows from the require- 
ment 

Performing the Lorentz boost to LF we obtain 

(1 - 3/3 2 ) T tt + 4(3T tz + (P 2 - 3) T zz = 0. (19) 

This is a quadratic equation for /3, which has two roots, 

/3i, 2 =o± Va 2 + b, (20) 

where 



2T* 



3T Z 



3T" - T z 



(21) 
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Figure 2. Left: the magnitude of transversal polarization II(i9; 7) due to up-scattering by a relativistic electron as a function of observing 
angle in the local co-moving frame. The case of locally isotropic ambient radiation field is shown for three different values of Lorentz 
factor 7. In the inset the emission diagram shows the corresponding lab-frame polarization. The lobes become gradually flattened toward 
the front direction of motion as 7 increases. One can see from eq. (11) that the graph of n($) is symmetrical with respect to 1? = 7r/2. 
This corresponds to a well-known fact that polarization is maximum for the radiation scattered perpendicularly to axis of symmetry in 
the CF. Right: Contours of II(i9, /3) are shown superposed on the density plot of 7(i9, 0). Levels of shading give the intensity (in arbitrary 
units) and illustrate the progressive beaming towards 1? = direction in the ultra-relativistic limit (in the lab frame). On the other hand, 
given a value of f3, the polarization degree II(i9; /3) as function of $ reaches maximum at a non-zero angle, always off axis (dashed line). 



Clearly, eq. (18) can be satisfied independently of the direc- 
tion of observation. For /3 — > /3i,2 the polarization changes 
from longitudinal to transversal. 

Secondly, we introduce the saturation velocity /3o 
(Sikora & Wilson 1981). As was shown by various au- 
thors under different approximations about the particle 
cross-section and the form of gravitational field (see e.g. 
Abramowicz et al. 1990; Vokrouhlicky & Karas 1991; Melia 
& Konigl 1989; Fukue & Hachiya 1999; Keane et al. 2001), 
the saturation velocity plays an important role in the 
dynamics of relativistic jets: particles moving at velocity 
smaller/greater than the saturation velocity gain/lose their 
momentum at the expense of the radiation field. In absence 
of other acceleration mechanisms and neglecting inertia of 
particles, the effect of radiation pressure eventually leads to 
— > Po as terminal speed of the particle motion. 

The saturation velocity is determined by the require- 
ment of the vanishing radiation flux in CF, i.e. 



0. 



This gives another quadratic equation, 



(1+/T)T 12 -/3(T" + T Z 
with the solution 

i-VT^2 

PO — 



0. 



(22) 



(23) 



(24) 



where a = 2T tz /(T u +T tz ). We ignore the second solution, 
as it has no physical meaning. 

As an example let us assume the incident radiation field 
to be strictly isotropic in the laboratory frame, i.e. 



T aP = diag (£ , \£ , \£ , \£ ) 



(25) 



with £ = T u being energy-density of radiation. Evaluating 
the stress-energy tensor in CF we find n m = -P 2 . Substi- 
tuting into the equation (13) we obtain polarization degree 



p 2 sin 2 d 



1+P 2 cos 2 § ' 
Lorentz transformation to LF gives 



n(tf,/3) 



(3 2 sin 2 ■& 



(2 7 2 - 1) (1 - P cos -&Y - P 2 sin 2 ■& 



(26) 



(27) 



Since II m ^ 0, the scattered radiation is polarized trans- 
versely. The critical velocities are Po — Pi — 02 = in this 
case. 

Figure 2 shows the dependence of II on observing angle 
according to equation (26). It is worth noticing that we ex- 
plore the frequency-integrated model, because this assump- 
tion is adequate for the purpose of clarification of the role 
gravitational lensing (discussed in the next section). The 
same dependency of Stokes parameters on the scatterer ve- 
locity and observer viewing angle is obtained in frequency- 
dependent calculation with the spectral index of the incident 
radiation equal to —1 (the case adopted originally by Begel- 
man & Sikora 1987). Our results are consistent with Lazzati 
et al. (2004) provided that appropriate averaging over en- 
ergy is adopted. It can be seen (in the left panel of Fig. 2) 
that the resulting curves closely resemble the numerical re- 
sult of Lazzati et al. (2004; cp. their Fig. 1). In particular, 
the curves are identical for 7 > 1 and they approach the 
ultra-relativistic limit II = (1 — cos 2 $)/(l + cos 2 (origi- 
nally examined again by Begelman & Sikora 1987, and in- 
voked more recently e.g. by Shaviv & Dar 1995). This limit 
corresponds to the case of a head-on collision, when all pho- 
tons are impinging at incident angles 7?i — * n because of 
aberration in CF. 

For moderate Lorentz factors there is some difference 
between our profile of LT(i9) and the corresponding numer- 
ical values plotted in Lazzati et al. (2004). For example, 
checking the 7 = 2 curve, we notice that relative difference 
amounts to roughly 13%. This apparent discrepancy is ex- 
plained by realizing that our eq. (26) has been derived in 
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terms of bolometric intensities, whereas Lazzati et al. em- 
ploy specific (frequency-dependent) quantities. By integrat- 
ing their Stokes parameters over frequency we recover pre- 
cisely the value predicted by eq. (26). See Horak (2005) for 
detailed comparisons. 



3 POLARIZATION OF LIGHT SCATTERED 
NEAR A COMPACT STAR 

3.1 The gravitational and radiation fields 

The gravitational field of a spherically symmetric star is 
described by Schwarzschild metric (Chandrasekhar 1992), 



ds 2 = -c 2 £di 2 + C 1 dr 2 +r 2 dO 2 , 



(28) 



where dfi is the angular part of a spherically symmetric 
line element, £(r) = 1 — Rs/r is the redshift function 
in terms of Schwarzschild radius, Rs = 2GMc~ 2 =2.95 x 
10 5 (M/Mq) cm, and M is mass of the star. Four- vectors 
and four-tensors will be expressed with respect to a local or- 

with non-vanishing 



thonormal tetrad, (e ( *\ e (r) , e (e) , e w 



-1/2 



and 



components e£*' = c£ 1//2 , e$. r ' 
e^f' = rsin#. Tetrad components of four- vectors are de- 
noted by bracketed indices and are raised and lowered using 
the Minkowski metric. 

Primary photons are emitted from a star and form the 
ambient radiation field acting on the particle. The star of 
radius R* and compactness R*/Rs appears to a static ob- 
server, located at radius r, as a bright disc of angular radius 
a* = a(r), 



sin a(r) = 



R e(r) 



1/2 



r ZiR) 1 / 2 ' 



(29) 



where R = max{|i?s, R*}- Because of light bending the 
solid angle subtended by a compact star on the sky is larger 
than the Euclidean (flat space) estimate. Formula (29) was 
originally discussed by Synge (1967) and the manifestation 
of the gravitational self-lense effect was examined by Win- 
terberg & Phillips (1973). 

In case of very high compactness (when R+ < |i?s) 
the rim of the image is formed by photons encircling the 
perimeter of the star more than once. In spite of compli- 
cated trajectories of photons, to the observer the surface 
appears radiating with intensity (we neglect limb darkening 
for simplicity) 



/(r) = 



h(R*). 



(30) 



Let us take the previous example from eq. (25), but assume 
that the source of primary photons occupies only a fraction 
of the local sky of the scattering particle. Gradual dilution of 
the source radiation with distance is described by function 
a(r). The limit of a — > n corresponds to strictly isotropic 
radiation arriving from all directions, whereas for a — > wc 
obtain the case of a point-like source. We thus recognize the 
results of subsection 2.3. 

The stress-energy tensor of the stellar radiation field 
has three independent components, namely, the energy den- 
sity, the energy flux, and the radial stress. These are given, 



2tt 



respectively, by 

ee Tl tm = ^-/(l- cos a), 
= cTi t)w =ttJ sin 2 a, 
V, ee TWW = ^7(l-cos 3 a). 



(31) 
(32) 
(33) 



There are two other non-zero components, — T+^^\ 

which can be computed from the condition T+ a — 0. Magni- 
tude of the star radiation is characterized by total luminos- 
ity, L-,, = AitRlT^R*). 

Finally we include another, isotropic component of the 
radiation field with intensity 7i so in addition to stellar light. 
The corresponding stress-energy tensor is entirely deter- 
mined by energy density £ lso = 47rc~ 1 7i so . The stress- 
energy tensor of the total radiation field is the sum T al3 — 
rpap _|_ Combining the two contributions allows us to 

model different situations according to their relative magni- 
tude and the motion of the scattering medium. 



3.2 Polarization of scattered light 

We first assume velocity of the scatterer /3(r) and compute 
the resulting polarization. Figure 3 shows the effect of van- 
ishing and changing polarization which occurs at particular 
value of [3. We compare two situations: the case of purely 
stellar component of the primary irradiation, as described in 
the previous paragraph (in the left panel), versus the case 
of a sum of the isotropic component and radiation coming 
from the star surface (on the right). The latter configuration 
represents an anisotropic irradiation of an electron; it can be 
parametrized by the mutual ratio of redshifted radiation in- 
tensities received at a distant observer, i.e. 



Ai ee — y^. 

/* 



(34) 



This can be considered as a toy-model of inverse Compton 
up-scatter in an illuminated jet where intensity of ambient 
light is not directly connected with the intensity of the cen- 
tral source. Ai = const is a free parameter of the model; given 
a value, the degree of anisotropy depends on the distance 
from the star in R§. It is worth noticing that light-bending 
is taken into account in this calculation automatically, in- 
cluding all higher-order images encircling the star. 

Polarization is non-zero provided that particle velocity 
is not equal to Pi,2(r) and, indeed, II can reach large values. 
This is shown in figure 4, where we plot the extremal value 
of function II m (/3, £) in the plane of particle velocity versus 
distance. |II m (/3, ()\ is equal to extremes of the polarization 
degree measured along a suitably chosen observing angle 
The curve of zero polarization is also plotted and we notice 
that it is independent of In this figure the primary un- 
polarized light was assumed to be a mixture of stellar and 
ambient contributions (the latter component was assumed 
to be distributed isotropically in the lab frame). The satu- 
ration curve A)(C) is shown and it is worth noticing that, 
for some values of the model parameters, /3o(C) crosses the 
contour of II = 0. Therefore, a hypothetical particle mov- 
ing or oscillating along the saturation curve would exhibit 
polarization that swings its direction by right angle. 
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Figure 3. Left: The case of incident radiation originating from an isotropic source of angular radius a; see eqs. (31)-(33). Two branches 
of critical velocity are shown, f3\(a) and /32(a), at which the total polarization of scattered light vanishes independent of observing 
direction. The saturation curve (3o(ct) is also plotted assuming that the radiation drag dominates the particle dynamics. Right: the 
same as on the left but for a mixture of two components of the incident radiation, i.e. the ambient isotropic (a = n) source plus stellar 
(non-isotropic) contribution according to eq. (34) with A; = 0.001. In both panels, the regions of longitudinal and transversal polarization 
are distinguished by shading. 



One can modify the previous example by considering 
a constant ratio of energy density, i.e. by replacing Ai with 
another parameter, 



Si* 



(35) 



This definition captures better the case when the ambient 
light originates from scattering of the central component 
(perhaps on clumps being accreted onto the star), so that 
both contributions are linked to each other and their en- 
ergy density decreases at identical rate with the distance. 
We again constructed graphs of II m (/3, C) an d found a sim- 
ilar structure of contours at small radii as those shown in 
Fig. 4, including the double-valued function /3i,2(C)- How- 
ever, the saturation velocity /3n(C) does not fall to zero at 
r 00, instead, it generally reaches substantially higher val- 
ues. Moreover, the critical point (where the separatrix curve 
self-crosses) is lost, as well as the whole structure towards 
right of the critical point. 

Polarization of scattered light obviously depends on 
scatterer motion and these can be calculated consistently. 
In the next section we finally determine velocity /3(r) along 
the particle trajectories, for which luminosity of the star and 
its compactness are parameters. 



3.3 Polarization along the particle trajectory 

Four-velocity it of a particle can be found by solving the 
equation of motion in the form (Abramowicz et al. 1990), 



VP a 
U 



(36) 



where m is the particle rest mass and K = + c~ 2 u"u li is 
a projection tensor. Left-hand side of eq. (36) includes the 
effect of gravity (V denotes covariant differentiation with 
respect to curved spacetime geometry) and the right-hand 
side provides the effect of radiation drag - accelerating or 
decelerating particles with respect to free-fall motion. 



Non-zero components of four-velocity are vS^ — cy and 
cyfi in the local tetrad of a static observer. 2 Equa- 



,M 



tion (36) takes the form 



dip, 



e /2 k 



(r) 



c 2 R s 
2r 2 



(37) 
(38) 



where dot denotes derivative with respect to coordinate time 
t and F\j is the radial component of the radiation force 



?(r) 



cr T 7 



[(1+/3 2 



T 



MM 



PIT 



n(*)(t) 



+ T 



MM 



(39) 



The effect of radiation on the motion is expressed by the first 
term in the parentheses on the right-hand side of eq. (37), 
whereas the other term can be considered as the contribution 
of gravity. Hence, the particle dynamics depends on relative 
strength of the radiation and gravitational fields. Because of 
redshift factor near a compact star these two influences do 
not obey the same simple Newtonian law, and so a rich set 
of possible results emerges. These can be parametrized by 
Eddington luminosity Le, which follows from the condition 
of zero acceleration for matter hovering at radius r — R*. 
The radiation force becomes F^ d — (cr T L*)/(47ri? 2 c) and 
the equation (37) with (3 = gives 



Le = 



2nmc 3 Rs 



a T £(R i 



U/2 ' 



(40) 



2 In this paper we consider purely radial motion. The same for- 
malism can be readily applied to more complicated motion of the 
scatterer, although it will then hardly be possible to solve both 
the motion and the resulting polarization analytically. Notice that 
the case of clumps orbiting in the plane of a black hole accretion 
disc was discussed by Pineault (1977), Connors et al. (1980) and 
Bao et al. (1997). These authors also pointed out that effects of 
general relativity could be discovered by tracing time variable 
polarization. 
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Figure 4. Contours of extremal values of the polarization function IT m (/3, Q for photons Thomson-scattered on an electron, moving with 
a given velocity (3 through a mixture of stellar and ambient diffuse light. Each panel captures the whole range of radii from r = R* (£ = 0) 
to r — » oo (£ = 1). The three rows correspond to progressively increasing luminosity parameter (34): (i) Ai = (top); (ii) Ai = 0.001 
(middle); and (iii) Aj = 0.1 (bottom). The left column is for a highly compact star with i?* = 1.01i?si the right column is for i?* = 10 3 iJs- 
Hence, the light-bending effects are significant on the left and negligible on the right. The curve of zero polarization IT = is plotted with 
a dashed line. Generally, if the star is sufficiently compact then the curve of zero polarization becomes double-valued; its two branches 
correspond to /3 = /3i,2(C) i n the previous figure. A separatrix is a particular contour that distinguishes regions of different topology in 
this graph. The saturation curve /3o(C) i s a l so plotted. 
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Figure 5. Top row: particle velocity /3(C) (thick curve) and the three critical velocities /3o, /3i and in the combined radiation and 
gravitational fields. Each trajectory starts from the star, R+ = 1.205-Rs (C = 0)j an d it proceeds towards infinity (£ = 1). Positive values 
of /9 correspond to an outflow, negative values are for an inflowing material. Three cases are shown for different values of dimensionless 
luminosity: T = 100 (left), T = 2 (middle), and T = 0.1 (right). Bottom row: the polarization magnitude 11(C) along the particle 
trajectories corresponding to the three above-given solutions. In each panel, the curves are labelled with the observing angle All curves 
have a common zero point where they cross each other. The sign of II distinguishes here the case of transversal polarization from the 
longitudinal one. 



We note, that the acceleration depends on the radial dis- 
tance from the center as well as on particle velocity. The 
relative importance of radiation and gravity is characterised 
by dimensionless factor 

?=T- ( 41 ) 

The radiation term in acceleration is regulated by the in- 
terplay of relativistic aberration and the Doppler boosting, 
which tend to establish the saturation velocity. At this point 
further acceleration vanishes, i.e. $o{r) = 0. Considering 
only radiation from the star and expressing the explicit form 
of the stress-energy tensor, eqs. (37)-(38) reduce to eq. (2.3) 
of Abramowicz et al. (1990). With gravitational attraction 
of the centre taken into account, a possibility occurs of an 
equilibrium point £ oq = f(r = R cq ), where a particle can 
reside. By setting (3 = and = in eqs. (37)-(39), one 
can find that the equilibrium radius ranges from i? cq (r)— >R+ 
(i.e. £ cq — > 0) for T — > 1 up to R eq (T) — > oo (i.e. C cq — > 1) 

for r -> Vs. 

Equations (37)-(38) allow for a finite set of topologically 
different solutions. These can be classified into different cat- 
egories (Abramowicz et al. 1990; see also Keane et al. 2001) 
according to the behaviour of saturation curves in (j3, em- 
plane. Notice that we already examined one of these solu- 
tions, i.e. the saturation curve (24) for very high luminosity 
of the star and negligible mass of the particle, i.e., F — > oo. 



The motion is then governed solely by radiation drag. We se- 
lect this condition because it is particularly relevant for the 
discussion of the resulting polarization of scattered light. Its 
role can be inferred also from Fig. 4, where the /3o-curve for 
this case passes through the critical point of contour lines of 
II(/3, £). The limit of V — > oo is an extreme case. Different 
profile Aj(C) applies to moderate values of the luminosity 
parameter, T < oo, when particles do not strictly maintain 
the saturation velocity because of inertial effects acting on 
them. 

Different categories of the particle motion then provide 
a natural framework also for the discussion of the resulting 
polarization. Three most important cases are recorded in 
figure 5. In this example only the stellar radiation is taken 
into account, whereas the component 7i so = for simplic- 
ity. The cases shown here correspond to the situation when 
(i) the radiation field dominates over gravity and the elec- 
tron is therefore pushed away to an infinite radius (see the 
left panel); (ii) a moderate value of the luminosity allows 
the scattering electrons to reach an equilibrium position at 
(" = C cq = 0.62 (middle); (iii) the luminosity is very small 
and the particles are almost free-falling in the gravitational 
field (right). Particles start from £ = and they quickly 
adhere to the saturation curve f3o(^;T), provided that ra- 
diation is dynamically important, i.e. in cases (i) and (ii). 
This occurs independently of initial velocity; then the mo- 
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tion follows a curve adjacent to but slightly different from 
the saturation curve. On the other hand, in the case (iii) the 
gravitation governs the motion; the trajectory /3(C) is only 
slightly asymmetric with respect to j3 = line by the weak 
influence of radiation. 

By coupling the equations of particle motion with the 
polarization equations of sec. 3.2 we obtain Stokes parame- 
ters of scattered light along each particle trajectory. Bottom 
panels of Fig. 5 show the resulting magnitude of polariza- 
tion. Notice how it crosses zero level at certain distance of 
the scatterer from the stellar surface. At this point polar- 
ization changes direction from transverse to longitudinal. 
Points of intersection of curves /3i,2(C) with the particle mo- 
tion /3(C) determine the radial location of the point of van- 
ishing polarization (indicated by dotted vertical lines in the 
plot). 

We will now assume that a small cloudlet is formed by a 
group of electrons. We can distinguish three cases, depend- 
ing on the bulk velocity of the cloudlet. These are discussed 
in subsections 4.1-4.3 below. The predicted time dependence 
offers a way to test the model. 



4 POLARIZATION FROM A CLOUDLET 

4.1 The case of fast ejection (7 3> 1) 

Let us denote R c \ >C Rs radius of the cloud and tp^R c \/z 
its angular radius as seen from the center. We assume that 
the cloud has small optical depth, t c i < 1, and it is ejected 
along z-axis, with bulk velocity f3(z) >• directed approxi- 
mately toward an observer (inclination angle of the observer 
is denoted O ). Clearly, 7 ^> 1 implies the scattered pho- 
tons are boosted in the direction of motion. Although light 
bending increases the apparent size of the star on the par- 
ticle local sky, general relativity effects are quite negligible 
on scattered photons moving straight away from the center. 
Only few photons are scattered backwards, and therefore 
the direct image greatly dominates the signal received by an 
observer. 

The measured radiation flux Stot has two components: 
the primary unpolarized flux 5 1 * and the flux of partially 
polarized radiation scattered in the cloud S c i. Their ratio 
can be given in terms of redshifted intensities and I c \ of 
the star and of the cloud, and by the ratio of solid angles 
occupied by the cloud and by the star on observer sky. This 
provides us with the estimation of the expected fractional 
polarization of the total signal. The ratio of fluxes is 



(44) 



s7 



Rci 

~r7 



id 



(42) 



Intensity arriving from the cloud I c \ = £ 2 (r)/ sc , where I BC {r) 
is the locally emitted intensity, as given by eq. (17). The 
scattered component is polarized with the magnitude II c i, 
as derived in the previous section. The total flux and the 
total polarization are 



Stot = (1 + s)S* 



lit 



slid 
1 + s' 



(43) 



Substituting equation (17) into (42) one can verify that 
the flux ratio s does not depend on the star intensity /*, and 
hence 



where k = r c i(i? c i/7?*) 2 depends on the size and density of 
the cloud, and / includes the geometry of the radiation field 
and the beaming/aberration effects arising from the cloud 
motion. We consider situations when k is small; then the two 
contributions to the radiation intercepted by the observer 
become comparable only in case of strong beaming, which 
leads to / 3> 1 for small observing angles. 

The moment of observation, i.e. the arrival time of pho- 
tons iobs( r — * 00 )i is related to the moment of emission t 
by 



t ohs ~ t — ^ — 2 j cos 9, 



(45) 



where we set t(zo) = for the initial time and t ohB — for 
the moment when the signal arrives to the observer. Notice 
that this estimate is sufficient for direct image photons dis- 
cussed in this subsection, but it would not be appropriate 
for higher-order image photons in the next subsection (in 
Schwarzschild geometry one can express time of arrival in 
terms of elliptic integrals, proceeding in the same analytical 
manner as above in the calculation of the ray trajectory; see 
also Bozza & Mancini 2004; Cadez & Kostic 2005). 

The temporal behaviour is shown in figure 6. We as- 
sume that the cloud has been pre-accelerated to large initial 
speed (3{t = 0) near the star surface. The graph captures 
the subsequent phase of gravitational and radiative decel- 
eration. The scattered light contributes significantly to the 
total signal only for a short initial phase (a peak occurs in 
the graph). The local maxima of the radiation flux (at t = ti) 
and of polarization (at t = tp) can be understood in terms 
of beaming: most of radiation from the cloud is emitted in 
a cone with the opening angle ~ I/7 about the direction 
of motion. For small viewing angle (0 o < 13 deg) the ob- 
server was initially located outside this cone but, as time 
goes, the electron decelerates, the cone opens up and the 
observer intercepts more radiation. The maximum observed 
polarization occurs with a certain delay tp — 1\ (proportional 
to M) after the peak of radiation flux. Subsequent decay of 
the signal is connected with a diminished scattering power 
of the cloud and the overall dilution of the radiation field. 
The observed polarization and the flux are lagged with each 
other and sensitive to the angle of observation. This be- 
haviour is clearly seen also in figure 7, where we assumed 
several different viewing angles. 

We selected large initial velocity in this example, oth- 
erwise the effects of aberration and the Doppler boosting 
would be less prominent, time-scales longer, and the effect 
of fractional polarization crossing zero point would disap- 
pear. The time span of this plot can be scaled according to 
the light-crossing time in physical units, 



1.5 ^ = 1.5 x HP 1 



M 



10M Q 



[secj, 



(46) 



i.e. proportionally to the central mass. The polarization 
magnitude is correlated with the intensity (this correlation 
was already noticed for the isotropic radiation in the right 
panel of Fig. 2). 
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Figure 6. Left: the variation of the total normalized radiation flux Stot(i) and the corresponding degree of polarization ITtot(t). The 
case of ultra-relativistic ejection starting with initial 7(i?*) = 10. Right: The corresponding velocity profile /3(C) and critical velocities 
fil, 2(C) are shown. Parameters of the plot are 8 = 17 deg, = 1.2i?g. Polarization vanishes at the point A when /3(C) = Al(0; the 
corresponding time is t = 4a in the left panel. Polarization vanishes once again at a later time, when /3(C) = /32(C)- 
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Figure 7. The relation between the normalized radiation flux and 
the total polarization. The magnitude of polarization Iltot reaches 
up to ~ 65% for suitable view angles. Values of o (in degrees) 
are given with the curves. Other parameters are the same as in 
previous figure. 



4.2 A cloudlet at rest (7 = 1) 

An interplay between gravity and the ambient radiation 
stalls the bulk motion, /3(t) — > 0. Scattered photons are 
then no longer boosted in the outward direction, and so the 
higher order (highly bent) rays can provide a non-negligible 
contribution to the observed light after encircling the star. 
This of course requires large compactness; we set R+ = |i?s 
hereafter. Again we assume an observer near z-axis and a 
cloudlet with a small size, ip <C 1. Unlike a more traditional 
application of the lense geometry, the cloudlet is placed at 
an arbitrary finite distance z = z(t) above the star and the 
deflection angle does not have to be small. 

Let us consider rays making a single round (by the an- 
gle O = 2-7T ± 6* ), with a radial turning point at pericenter 
r = r p tightly above the photon circular orbit. As mentioned 



in subsect. 3.3, the equilibrium radius depends on the star 
luminosity. Once the cloudlet settles at the equilibrium point 
r = -Rcq, scattered photons are no more boosted to high en- 
ergy and the collimation effect disappears. In this situation 
relatively more light is backscattered in the direction toward 
the photon orbit. Part of these photons form a retrolensed 
image (Holz & Wheeler 2002), which may also reach the ob- 
server. We thus now calculate (de)magnification of light also 
for the two first-order images, which give the most signifi- 
cant contribution and may influence the net polarization at 
infinity. To this aim we need to consider rays starting near 
above the star, passing through pericenter and eventually 
escaping to infinity (the retrolensing geometry; see Ohanian 
1987; Virbhadra & Ellis 2000; Bozza 2002, and references 
cited therein). 

Two arcs are formed which merge together in the Ein- 
stein ring (with radius just above b c — | \/3Rs) if the 
observer is aligned with the source. By integrating null 
geodesies, expanding the elliptic integrals near the pericen- 



ter r p ~ |^s, assuming the deflection angle close to O ~ 2ir, 



and keeping only the leading terms we obtain the desired 
width of the two retrolensing images, 

2 . 2 A 1/2 



5b(<j)) = 2 Sb (V> " 80 sin 

where 8b = K(z) e - 2jr y \4 
and 

3^3 %/E 



s3 



K = 



(47) 

^ arcsin(-!/j/6> ), ■ip — Min{0 o , tp}, 



3u 



, ^ _ Rs 
u(z) = — . 
z 



(48) 



2 ^+1 \/3 + \/l + 3u' 
We remark that O ~ 2-k was assumed for simplicity only. 
The case of arbitrary Q can be treated in similar way. 

Time dependency of the arcs is caused by the scatterer 
motion, z — z(t). We integrate over the cross-section of the 
arc images to derive their total luminosity and polarization 
at each time moment. Higher-order images suffer from the 
de-magnifying influence of the light bending, which reduces 
their luminosity, unless a special geometrical alignment of 
the source and the observer occurs and favours the oppo- 
site effect of a caustic. This can be quantified by the gain 
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Figure 8. Left: the gain factor M(&, f) according to the approximation formula (51). Two curves are shown as a function of the source 
distance (, for = 358° (solid line) and © = 362° (dashed line). For comparison, exact (numerically computed) values are also plotted 
with triangles and circles. Right: normalized intensity of the direct image (thin solid line) and the two retrolensed images (dashed and 



dotted lines) of light scattered from a particle residing in the equilibrium point 



i? cq , as a function of the Eddington parameter. 



Notice that the contribution of the two retrolensed images is almost identical and it amounts up to 
line). Inclination O = 2° in both panels. 



20% of the total signal (thick solid 



factor, A4, which determines the ratio of fluxes received in 
retrolensed/direct images. The problem translates to eval- 
uating the ratio of solid angles, Ai = jfjS where indices 
"i/o" refer to the angular size of the source with/without 
taking the light bending into account. In the case of a small 
(but finite) size cloudlet, we find 

M(z) = 6 \/3 A"(z) e" 2 "?/.- 1 A(0 O /</>), (49) 
where the term 

A(k) = -E(<f>,k), $ = arcsin [Min {AT 1 , ljl , (50) 

7T 

arises from the integration over the Einstein arcs. 
For tjj <C 9 the gain function is 

Formula (51) reduces to eq. (21) of Ohanian (1987) for z — * 
oo, £(z) — ♦ 1. In our situation, eq. (51) requires that the 
cloudlet is sufficiently small in size and its motion is directed 
somewhat sideways with respect to the observer view angle. 
Figure 8 compares M(C, 0) with the corresponding result of 
a numerical integration, showing that the approximation is 
sufficiently accurate for our purposes. 

Adding the contributions from different parts of the 
source has a depolarizing effect on the final signal, which 
we illustrate in figure 9. For the total magnitude of polar- 
ization of the retrolensing images we find 

n ret = n(iV)K<VV>), (52) 

where fl(i9 p h) is the polarization magnitude of light scat- 
tered in the direction towards the photon circular orbit and 

p(k) = , 2 , / cos 2<j>Jl- k 2 sin 2 (j> dcj>. (53) 
TA(fc) Jo v 

Functions A and $ were defined in eq. (50). Function p de- 
termines the shape of retrolensing images in the observer 
plane; see figure 10 (we resolve a narrow trace of these arcs 



on observer sky by enlarging their separation 5b = b — b c 
from the critical radius b c ). Polarization vectors of all three 
images have the same orientation, but they experience dif- 
ferent time delays and lensing along each trajectory. The 
contribution of the retrolensed images is now evident, and 
quite significant. Notice that the angle i? p h is the apparent 
angular size of the photon orbit as seen on the local sky of 
the cloudlet. It enters in eq. (52) because higher-order im- 
ages are formed almost exclusively by light scattered on the 
photon circular orbit. In our case, i? p h = a*(z). The polar- 
ization magnitude drops sharply if the observer inclination 
is less than the angular size of the cloudlet. 

4.3 Comparison between an outflow and an inflow 

Now we consider an intermediate situation with moderate 
velocity of the bulk motion (both an outflow or an inflow, 
i.e. /3 ^ 0). For a moderate outflow velocity, the result is 
shown in figure 11. We consider a particle on the decel- 
erating branch of the trajectory in a weak radiation field, 
r — > 0, which eventually reaches the turning point (and 
starts falling afterwards). The two retrolensed images con- 
tribute about 10% of the scattered flux at maximum and the 
trajectory crosses j3 = curve, where the polarization 

vector swings its direction. The outcome is quite different 
for matter infalling onto the star, because scattered photons 
are boosted in the downward direction and a considerable 
amount of light is then directed on the photon orbit. As a 
result, the retrolensed images are more pronounced and they 
cause a brief flash of light. The resulting signal is shown in 
figure 12. 

The effect of retrolensed images is clearly visible in the 
polarization curve; see figure 13. The case shown in the left 
panel exhibits a brief drop of the polarization magnitude 
when velocity crosses /3 = /32(C) (the direct image arrives 
at t ~ 1 in dimensionless units). At this moment the polar- 
ization changes its orientation between transversal and the 
longitudinal one. Then the signal restores back to a non-zero 



On the role of strong gravity in polarization 13 





-10 o 10 
e [cleg] 



40 



Figure 9. Polarization magnitude from scattering on a particle at rest at z = i? cq . This represents a cloudlet of angular radius ip on the 
local sky of the star. Left: the observed polarization magnitude as a function of Eddington parameter T. Notice that R eq is a function 
of r, and so the graph covers the whole range of radii from the star surface to infinity. The observer inclination is O = 2°. Right: the 
corresponding polarization magnitude for different inclinations and constant T = 1.6 (in case of precise alignment, o = 0°, polarization 
vanishes because of symmetry). 
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Figure 10. Left: the form of Einstein arcs (a— c) and the ring (d) corresponding to the retrolensing images in polar coordinates (ft, rf>) in 
the observer plane. The source is supposed to be a circular target of angular radius ip = 2°, located on z-axis at distance z = 3Rs- The 
observer is at r — * oo and she has a small angular offset from the perfect alignment - (a) O = 10°, (b) O = 2°, (c) O = 1°, (d) O = 0°. 
Right: a contribution to the polarization produced by the Einstein arcs. A detail of the normalized magnitude p is plotted for small 
values of inclination; see eq. (52) for the definition of function p(0 o /ip). For large O the magnitude of polarization saturates at roughly 
constant level, equal to the polarization scattered in the direction toward the photon circular orbit. Two curves are parametrized by the 
angular size ip of the cloudlet, as indicated in the plot. 



value and the same behaviour repeats when the retrolensed 
image arrives after a certain delay (at t ~ 21). The case 
shown in the right panel exhibits a similar flash, also caused 
by the contribution of the retrolensed photons. However, 
now we observe a single fluctuation, which is actually an 
increase of the polarization magnitude; this is because the 
case shown here corresponds to transversal polarization dur- 
ing the whole observation and the trajectory does not cross 
any of the critical curves = /3i,2. 



5 CONCLUSIONS 

Our calculation here is self-consistent in the sense that the 
motion of the blob and of photons, and the resulting polar- 
ization are mutually connected. We concentrated ourselves 
on gravitational effects and neglected other intervening pro- 
cesses, first of all the effect of magnetic fields to which po- 
larization is sensitive (see e.g. Agol & Blaes 1996). This al- 
lowed us to compare polarization magnitudes of direct and 
retrolensed images, which could point to the presence of a 
highly compact body. We have noticed the mutual delay be- 
tween the signal formed by photons of different order. The 
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Figure 11. An outward-moving particle decelerating in the gravitational and weak radiation fields. Left: the trajectory in the (/3, £)- 
plane, i.e. dimensionless velocity versus distance (thick solid line). Velocity changes in the direction of an arrow. Right: the radiation 
flux of scattered light as a function of time for the direct and two retrolensed images. The initial condition is fi = 0.7 at r = 2Rg. The 
observer inclination is O = 5°. 
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Figure 12. The same as in Fig. 11, but for fast inward motion of the scatterer and strong radiation of the star (T = 10). We set j3 = —0.99, 
r = 20-Rs as an initial condition at t = 0. The signal of the higher-order image flashes for a brief moment around dimensionless time 
t ~ 60. The intense radiation of the star reverses the particle velocity to the outward motion at later stages. The observer inclination is 
0o = 2°. 



delay is characteristic to the effect and has a value propor- 
tional to the central mass. 

Polarimetric properties are susceptible to large changes 
depending on detailed physics and geometry of the source, 
and this is not only the advantage, which could help us 
to trace how different objects are functioning, but also a 
complication. In particular, the polarization is sensitive to 
the source orientation and its magnitude fluctuates from 
case to case. Our results are useful for testing more com- 
plex and astrophysically realistic models with up-scattering 
of soft photons in jets and fast flows around black holes. 
In the black hole case the primary photons would be pro- 
vided by an accretion disc rather than the star surface and, 
hence, one can no more take advantage of its spherical sym- 
metry, which helped us to simplify our calculations here. 
Nonetheless, the same formalism can be employed and sim- 
ilar features of time-dependent polarization lightcurves are 
expected; strong gravity plays a vital role again. Putting 



this in another way: provided a 'realistic' equation of state 
implies neutron star radii greater than the photon circu- 
lar orbit, detection of the signatures of retrolensing images, 
which we discussed above, would exclude a neutron star as 
a candidate on the central body. 
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Figure 13. The magnitude of polarization corresponding to the case shown in Fig. 11 (left panel) and in Fig. 12 (right panel). The 
retrolensed signals are delayed with respect to the direct signal. The delay has a value characteristic to light-travel time along the photon 
circular orbit (it scales proportionally to the mass of the central body). The total polarization is suppressed or enhanced depending on 
the mutual relation between the polarization of direct and retrolensed photons. 
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